Dark-bright gap solitons in coupled-mode one-dimensional 

saturable waveguide arrays 

Rong Dong, Christian E. Riiter, Detlef Kip 

Faculty of Electrical Engineering, Helmut Schmidt University, 22043 Hamburg, Germany 

Jesiis Cuevas 

Grupo de Fisica No Lineal, Departamento de Fisica Aplicada I, 
Escuela Politecnica Superior, C/ Virgen de Africa, 7, 41011 Sevilla, Spain 

Panayotis G. Kevrekidis 
Department of Mathematics and Statistics, 
University of Massachusetts, Amherst MA 01003-4515, USA 

Daohong Song, Jingjun Xu 
The Key Laboratory of Weak-Light Nonlinear Photonics, 
Ministry of Education and TEDA, Applied Physics School, 
Nankai University, Tianjin 300457, China 

Abstract 

In the present work, we consider the dynamics of dark solitons as one mode of a defocusing 
photorefractive lattice coupled with bright solitons as a second mode of the lattice. Our investiga- 
tion is motivated by an experiment which illustrates that such coupled states can exist with both 
components in the first gap of the linear band spectrum. This finding is further extended by the 
examination of different possibilities from a theoretical perspective, such as symbiotic ones where 
the bright component is supported by states of the dark component in the first or second gap, or 
non-symbiotic ones where the bright soliton is also a first-gap state coupled to a first or second gap 
state of the dark component. While the obtained states are generally unstable, these instabilities 
typically bear fairly small growth rates which enable their observation for experimentally relevant 
propagation distances. 



PACS numbers: 05.45.Yv, 42.65.Tg, 42.65.Jx, 42.65.Hw, 42.82.Et, 63.20.Pw 



I. INTRODUCTION 



The examination of the Hamiltonian continuum model with periodic potentials and its 
discrete analog of lattice dynamical systems has been a topic of increasing popularity over 
the past few years [Ij. This is mainly due to their wide applicability in diverse physical 
contexts including, but not limited to, the spatial dynamics of optical beams in coupled 
waveguide arrays [2], optically-induced photonic lattices in nonlinear optics [3J, temporal 
evolution of Bose-Einstein condensates (BECs) in optical lattices in soft-condensed matter 
physics [1], and the DNA double strand in biophysics [3]. 

A principal research theme in this direction is the study of existence and stability of 
coherent structures in these models and their feasibility in experiments. Several years ago, 
fabrication of nonlinear optical AlGaAs waveguide arrays [6J provided a first prototype 
where many initial investigations arose, such as discrete diffraction, Peierls barriers, diffrac- 
tion management ^j, and gap solitons [8]. So far numerous fundamental investigations have 
been pursued in waveguide arrays including modulational instability [U], four- wave-mixing 
effects arising from the coupling of multiple components [TU], as well as the study of interac- 
tions of solitary waves with surfaces [H] . Subsequently, the formation of optically-induced 
photonic lattices in photorefractive crystals became an ideal platform for the observation of 
various types of solitonic structures. The theoretical proposal |12] and rapid experimental 
realization of such (mainly 2D) lattices [ISIIII], enabled the observation of, among others, 
dipole [U], necklace [ISj, and rotary [T7] solitons as well as discrete [ISl US] and gap [20] 
vortices. Recently, waveguide arrays in lithium niobate (LiNbOs) crystal, which possess a 
self-defocusing nonlinearity, have found significant applications in the study of modulation 
instability [21], beam interactions [22j, dark discrete solitons [23], bright gap solitons [21], 
dark solitons in higher gaps [25], as well as Rabi oscillations [26j. 

Our goal in this work is to consider the case of vector solitons. Although they have been 
studied both in the focusing case of bright- vector solitons in strontium barium niobate [27] 
and the defocusing case of bright-gap- vector solitons in LiNbOs [28] , much less work has 
been done in multi-component settings. Instead of mixtures of two solitary waves of the 
same type as in the above cases, we aim to examine the mixture of a bright with a dark 
soliton in photorefractive defocusing waveguide arrays. Such dark-bright states were first 
created in the absence of lattices in photorefractive crystals over a decade ago [29J and their 
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interactions were partially monitored [30] . In the context of BECs such solitary waves were 
also predicted theoretically |3T] , and generalizations thereof were considered as well (such as 
e.g. the dark-dark-bright or bright-bright-dark spinor variants of [32j). However, it was only 
quite recently that such structures were experimentally observed [33ti36] . This has led to a 
renewed interest in this theme, by addressing the interactions of dark-bright solitons from 
an integrable theory [37] or numerical [38j perspective, as well as their higher-dimensional 
generalizations [39j. To the best of our knowledge, there is no earlier investigation of such 
states in models with a periodic potential except in the context of nonlinear dynamical 
lattices [iU] . 

Our motivation, presented in section II stems from an experiment in defocusing LiNbOs 
waveguide arrays where a dark soliton state in the first gap (we will refer to this type of state 
as "bubble" in what follows) is coupled to a bright soliton in the same gap. We will show 
that these two waveforms coexist as a solitonic entity. Also, we will present conditions under 
which such a molecule may break up in its constituents. This, in turn, motivates a more 
detailed theoretical study of the different types of dark-bright states that can exist in the 
system. Such coupled states will be identified between either a bubble (in the first gap) or a 
higher-gap (i.e., the second gap in this case) dark soliton in the one component with either 
a regular bright soliton or with a bright gap soliton. When a bubble or dark soliton couples 
to a regular bright one, we refer to these solitons as symbiotic because the bright component 
can not exist without the supporting dark component (due to the defocusing nature of the 
nonlinearity) . For the coupling with a bright-gap soliton, because both components can 
persist individually, we refer to these states as non-symbiotic. In section III, we set up the 
model problem and benchmark it against experimental data by identifying its linear band 
spectrum. In section IV, the numerical results for the above soliton families will be given. 
Finally, in section V we summarize our findings and present conclusions as well as some 
relevant directions for future study. 



II. EXPERIMENTAL MOTIVATION 



To experimentally investigate such molecular solitonic states of dark and bright solitons, 
we used a ID waveguide array (WA) fabricated on an iron-doped lithium niobate (LiNbOs) 
substrate by in-diffusion of titanium at high temperature. Arising from the bulk photovoltaic 
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effect, the substrate crystal displays a saturable type of defocusing nonlinearity [H]. The 
transverse direction z is parallel to the ferroelectric c-axis. The direction of light propagation 
is along the y-axis. The array investigated in the following experiments consists of 250 
channels and has a grating period A = 8.5 /im, which is the summation of the channel width 
of 5 yum and a spacing of 3.5 /xm between adjacent channels. One of the end facets of the 
waveguide array sample is polished to optical quality to allow for direct observation of the 
out-coupled light from the array with the help of a CCD camera. 

In our experimental setup, we employed the prism-coupler scheme, with which we can 
selectively excite different Bloch modes in any desired band. Furthermore, with this method 
we can determine accurately the band structure of the waveguide array [12]. The experi- 
mental layout is sketched in Fig. [ij First, the input light with a wavelength of 532 nm from 
a frequency-doubled Nd:YV04 laser is expanded by a beam expander into a plane wave and 
then split into two separate beams. One beam propagates through a phase mask covering 
half of the beam along the transverse direction z. As a consequence, the covered half of 
the input beam experiences an additional vr phase shift, thus a dark notch is generated 
at the center of the intensity profile. Another beam is modulated by an oscillating mirror 
driven by a function generator. With applied external modulation, this beam is mutually 
incoherent with respect to the other beam. With the combination of two cylindrical lenses 
LI and L2, the beam passing the phase mask is then imaged onto the waveguide. Here the 
focal lengths of the two lenses are chosen in order to generate an ideal width of the dark 
notch covering about two channels, which is the input light pattern for the excitation of a 
dark soliton. The other beam is focused meanwhile by lens L2 with a diameter of roughly 
10 yum and serves as the excitation light for the bright soliton. Both beams are coupled into 
the waveguide array and co-propagate until they reach the end-facet of the sample. With a 
high resolution CCD camera, in combination with a 20x microscopic objective lens, we can 
monitor around 25 channels of the intensity distribution on the end-facet. With this setup, 
it is possible to adjust the input light distribution for both, the bright soliton and the dark 
soliton separately, for example, the relative locations of the two solitons on the waveguide 
array as well as different excitation angles for modes originating from different bands. 

In the experiment, a bright gap soliton was excited from the first and a dark soliton 
from the second band (a "bubble" according to our notation above), both at the edge of 
the Brillouin zone. The centers of both solitons were carefully adjusted to overlap on the 
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FIG. 1: Schematic experimental setup: HW, half-wave plate; P, polarizer; BE, beam expander; M's, 
mirrors; BS's, beam splitters; OM, oscillating mirror; FG, function generator; PM, phase mask; 
S, screen; LI and L2, cylindrical lenses; MO's, microscopic objectives; CCD's, CCD cameras; PD, 
photodiode; WA, waveguide array. 



same waveguide channel. We first checked under low optical power (less than 2nW per 
channel) the linear diffraction behavior of both the dark component [Fig. |2]^a), top row] 
and bright component [Fig. [2]^a), bottom row]. Then, by blocking one of the input beams, 
we formed individual gap solitons (either dark or bright) by increasing the optical power to 
appropriately high values [Fig. [2|^b)]. In all nonlinear experiments, the dark soliton from 
the second band was formed under 150 nW optical power per channel. In order to analyze 
the existence interval of the bubble-bright composite solitons, the input light power of the 
bright soliton was varied, resulting in different power ratios of dark and bright components. 
At first, a bright soliton was formed at 200 nW per channel, yielding a power ratio to the 
bright and dark solitons of 4:3. In this case, we observe a robust co-existence of the two 
components at the output facet, as shown in Fig. ^c). However, when we excite the bright 
soliton at much higher power (400 nW per channel, resulting in a power ratio of 8:3), the 
propagation constant fib of the bright component in this scenario is further decreased below 
the existence threshold (see also the theoretical analysis below), while the propagation 
constant fid of the dark component's bubble state remains essentially unaffected. The 
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FIG. 2: Experimental results showing linear diffraction (a) and soliton formation of individual 
bright and dark components (b). When both beams are launched with a power ratio of 4:3, a 
robust bubble-bright soliton is formed (c) where both, bright and dark components are centered 
on the same channel. 

result of the experiment in this situation is a clear spatial shift of the bubble center by one 
waveguide channel [Fig. [3|^a)] due to the coupling with the bright soliton. This shift may be 
understood as the initial phase of a repulsive interaction of the two constituents, and thus 
suggests the non-existence (or strong instability) of bubble-bright solitons for these input 
conditions. After reaching the steady-state for the input power ratio 8:3, we blocked the 
input beam used for excitation of the bright soliton. Because the nonlinearity in lithium 
niobate is non-instantaneous, the negative defect formed by the bright beam is still present 
and is only slowly erased due to the photoconductivity generated by the remaining dark 
beam. As a consequence, in the plane (see the left panel in Fig. [6]) we now move 

upwards (i.e., fib increases) on a vertical line, reaching back the existence regime of robust 
bubble solitary waves. We thus observe a reversible effect, presented in Fig. |3](b): namely, 
after the bright soliton is blocked, the dark soliton is restored to its original location. This 
restoration proves directly the repulsive influence from the dominant bright soliton. When 
the bright component is switched on again in Fig. ^c), once again the strong repulsion 
between bright and dark components forces the dark soliton to be shifted by one channel. 
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FIG. 3: When the input power ratio is increased (dominating bright component) a shift of the 
dark sohton center is experimentally observed (a). A restoration of the center position of the 
bubble state appears when the bright component is blocked (b). This process is reversible in the 
experiment, and the bubble is forced to shift again from the center position when the bright beam 
is switched on again (c). 

III. MODEL SETUP 

In what follows, we will consider composite solitons with a dark (or bubble) wave in one 
component coupled with a bright mode in the second component in the context of TE-TE 
modes for the geometry of our waveguide array. We start by presenting the underlying 
model in the full dimensional form with the paraxial approximation, and then we discuss 
the non-dimensional variant of the model which will be used for our numerical computations. 

A. Dynamical equations 

The paraxial equations for coupled TE-TE modes of the two beams represented by Ed 
and Eb in what follows, are given by: 

^^.E. + l^e..E. + ^ (An(Z) + A„„,^^|i^) = . 

with An{Z) being the refractive index profile and the propagation direction denoted as 
the x-direction. One can find "stationary" solutions of this system by defining: 



E,{X, Z) = e'^'^^'uiZ), E,iX, Z) = e*^''^i;(Z), (2) 

where Pdfi are the propagation constants in the X-direction and u{Z) and v{Z) the 
amphtude profiles of each TE mode, which, in turn, satisfy: 



-hv + + A + A„„,^^^) „ = , (3) 

The values used in the experiments are the following ones: 



Us = 2.2341, A = 532 nm, A = 8.5 /im, k = = 26.386 ^m-\ Auni = 2.5 x 10~^ 

A 

(4) 

(cf. also the discussion given in section II) where Ug is the refractive index of the LiNbOs 
substrate for extraordinary polarized light, A is the wavelength of the input light, A is the 
period of the waveguide array and Auni is the maximum refractive index change induced by 
the nonlinearity. 

The refractive index profile can be determined by adjusting the experimental Bloch bands 
showing the change of the effective refractive index Arics = fics — "^s, with = Pko and 
kf) = k/ug being the transverse wavevector in vacuum. The refractive index is then given 
by: 

An{Z) = Ano + AniV{Z) (5) 

with 

V{Z) = cos ( — j - 0.25 cos ( — ) (6) 

and 



Ario -ns = 27.567 x 10"^ Ani = 8.35 x 10"^ . (7) 

Figure |4] shows the correspondence between the experimentally observed Bloch bands [12] 
and the theoretically computed ones. Clearly, the above set of parameters offers a very good 
handle on the linear part of the problem. 
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B. Non-dimensional equations and parameters 



The non-dimensional version of the system of Eqs. ([T]) is given by 



id^u H — dzzU + [?7° + riV{z)]u + v 



1 + m2 + t;2 

id^v + -S^^f + \r]^ + r]V<yz)\v + v— — — — - 

2 1 + + f"' 



rM = 



f = 



(8) 



while the stationary states are solutions of: 



-IJ-dU + -d^zU + + r]V{z)]u + u— — — — -u = 
/ 1 + + 



t; = . 



(9) 



The non-dimensional parameters are related to the experimental ones by the following 
relations: 



u = ± 



a 



2 ' 



9 ' 



(10) 



(11) 



0°°°°°°°°°°°°°°°.-- ,, 




-Q.5 0.5 



FIG. 4: Bloch bands numerically (full lines) and experimentally (circles) determined. 



The parameter a has been introduced so that the non-dimensional values are of 0(1). 
Throughout the calculations, it has been fixed to a = 10. The sign of v indicates either self- 
focusing (positive) or self-defocusing (negative). Additionally, the nondimensional distances 
are given by 

z = aZ/A, X = = . (12) 

The refractive index profile and parameters are given now by 

V(.) = cos(—) -0.25 cos (—) (13) 

ri = 0.1880, r/° = 0.6207, u = ±0.0563 (14) 
and the change of the effective refractive index is 

for each (dark and bright) component. 



C. Stability equations 

Once stationary solutions of the boundary value problem (with periodic / anti-periodic 
boundary conditions, depending on the nature of the examined solution) of Eqs. (|9| are 
identified, their linear stability is considered by means of a Bogolyubov-de Gennes analysis. 
Namely, small perturbations [of order 0(5), with < 5 ^ 1] are introduced in the form 



Ediz,x) = e'^''^ [uoiz) + 6iPiz)e'^' + Q*(^)e-^-*^)] , 

Eb{z, x) = e^^"^' [vo{z) + 6{R{z)e'''' + S*{z)e-''^*')] , (16) 

and the ensuing linearized equation are then solved to 0(5), leading to the following 
eigenvalue problem: 

/ P{z) \ / Li L2 La L3 \ / Piz) \ 

Q{z) -L2 -Li -L3 -La Q{z) 

u = , (17) 

R{z) L3 L3 L4 L5 Riz) 

\ S{z) ) \ -L3 -L3 -L5 -L4 / V S{z) ) 
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for the eigenfrequency u and the associated eigenvector {P{z),Q{z), R{z), S^z))"^, where 
Lj, j = 1 . . .5 are the following operators: 



L3 



1 cf 

~f^d+-^ + [v' + vV{z)] + iy 



ul + vl 



+ 



Un 



i + ui + vi ' {1 + ui + viy 



Un 



UqVq 

,2\2 ' 



:i+ui 



7]V{Z)] + V 



+ 

\ + ul + vl 



'\ + uI + Vq 



2\2 



,2^2 ' 



(18) 



where it has been taken into account that uq{^z\ Vq(^z) G M. Once the stationary solutions 
are found to be linearly unstable (i.e., Im{u;} 7^ 0), then the dynamical manifestation of 
the corresponding instabilities is monitored through direct numerical simulations of Eq. (|8]). 
As we will see in the next section, all of the analyzed solutions are unstable, although their 
growth rates are so small that long propagation distances x are needed in order to observe 
the emergence of the pertinent instabilities. 



IV. NUMERICAL RESULTS 



We now present our results for the several types of coherent structures considered in our 
system in the self-defocusing setting (i.e. < 0). All of them are composed of a bright 
soliton in the 1st band gap. The dark structure can be of two types. It may be a bubble, 
located in the 1st band gap and arising from the top of the second Bloch band in which case 
the overall phase shift between the two endpoints of the domain is 0. Alternatively, it may 
be a (genuine) dark soliton, which emerges from the bottom of the second Bloch band, and, 
consequently, its propagation constant is found in the second band gap and it bears a phase 
shift of vr between the domain endpoints. 

We make one more terminological distinction between the different types of waveforms 
that can arise. In particular, the emerging bubble/dark-bright structures can either be 
symbiotic or not. In the first case, the bright soliton is unstaggered and emerges from the 
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top of the first band (zero mode). These modes are called symbiotic because an isolated 
bright component would not exist in this form for the relevant values of the propagation 
constant; it necessitates the formation of an effective potential by its dark (or bubble) 
counterpart in order to co-exist with it. In the second (non-symbiotic) case, the bright 
soliton is staggered and emerges from the bottom of the first band as a genuine gap soliton 
that would be sustained in the system even in the absence of the other component. 

These two distinctions (dark or bubble waves for the first component, symbiotic or non- 
symbiotic ones depending on the nature of the second component) give rise to four possibil- 
ities for the ensuing structures dubbed as follows: symbiotic / non-symbiotic bubble-bright 
soliton (SBBS / NSBBS) and symbiotic / non-symbiotic dark-bright soliton (SDBS / NS- 
DBS). Among the four, it is the NSBBS that was observed in our experimental motivation 
in section [Tlj Figure |5] shows prototype examples of the input field profiles for each of these 
four solutions. 

As mentioned above. Fig. |4] shows the position of the Bloch bands which are also relevant 
for the identification of the nonlinear localized modes that arise in the system. In particular, 
the first band is located in the interval fiu = 0.6755 < fi < 0.6833 = [3.00 x 10^^ < 
rics < 3.03 X 10"^], the second one is /i2d = 0.4614 < jj, < 0.5181 = /i2n [2.05 x 10"^ < Ucs < 
2.30 X 10-3] and the third one at /isd = 0.1806 < n < 0.3567 = fi3u [0.80 x lO^^ < < 
1.58 X 10-3]. 

Our numerical computations show that, in absence of coupling between the modes, 
the bright soliton can be identified in the first gap for fii^ + u < fi^ < fiia (i-e. fib G 
[0.6192,0.6755], and u here as well as below denotes an appropriate shift) whereas bubble- 
type solutions also exist for ^2u + < fid < fi2u (i-e. /id G [0.4618,0.5181]) and, in turn, 
the dark soliton can be identified for lower values of the propagation constant, namely for 
/^3u + ^ < fid < fi3u (i-e. /id G [0.3004,0.3567]). In the case of the two coupled beam 
components within the waveguide array, the existence interval is narrower. Furthermore, 
the existence range depends qualitatively on the symbiotic / non-symbiotic character of the 
soliton. More specifically, the accessible range of fid, for a given fif,, is always wider for sym- 
biotic solitons than for non-symbiotic ones. Additionally, the existence range of symbiotic 
solitons is limited from above by /ii^. Figure |6] depicts the existence range for dark-bright 
and bubble-bright symbiotic as well as non-symbiotic solitary waves. 

We have examined the linear stability of the obtained solutions, finding that the relevant 
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Z [mm] Z [mm] 

FIG. 5: Profile of the electric field for a NSBBS with fib = 0.64 [Anes,b = 2.84 x 10"^] and 
fid = 0.49 [Aneff,d = 2.18 x lO"^] (left panel in top hne). Profile of the electric field for a SBBS 
with Hb = 0.66 [Aneff,b = 2.93 x 10"^] and fid = 0.48 [Aneff,d = 2.13 x 10"^] (right panel in top 
line). Profile of the electric field for a NSDBS with fib = 0.65 [Ancg.b = 2.89 x 10"^] and fid = 0.33 
[Aricfr.d = 1.47 x 10^^] (left panel in bottom line). Profile of the electric field for a SDBS with 
fib = 0.67 [Aricff^b = 2.98 x 10"^] and fid = 0.32 [Aneff,d = 1.42 x lO^^] (right panel in bottom 
line). Blue dashed lines: input field of the bright component. Black solid lines: input field of the 
dark component. The red solid lines in each case illustrate a rescaled form of V{x) to indicate the 
location of the potential wells. 



waveforms are generically unstable in the spectral sense. I.e., we have identified an imaginary 
or complex eigenfrequency associated with the linearization spectrum around these profiles, 
however the growth rate is typically fairly small (< 10^'^ in non-dimensional units, i.e., 
< 0.05 mm~^ in dimensional units) and always less than 10~^ in non-dimensional units 
corresponding to 0.5 mm~^ in dimensional ones. Consequently, instabilities appear at a 
sufficiently large propagation distance X (inversely proportional to the above growth rate). 
Figure [7] shows the growth rate dependence with Anefr,d and Anes,h for the four analyzed 
structures; notice the colorbar on the right indicating the magnitude of the respective growth 
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FIG. 6: Existence range for bubble-bright (left) and dark-bright (right) solitons. Relevant end- 
points of the linear spectrum (and cutoff points below which we were unable to continue the 
solution) are denoted by corresponding horizontal or vertical dashed lines. 





FIG. 7: Growth rates for bubble-bright (top) and dark-bright (bottom) solitons. Left (resp. right) 
panels correspond to symbiotic (resp. non-symbiotic) structures. 
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rates. In order to test the effect of instabilities, a random perturbation of magnitude ~ 10~^ 
is introduced to the input field profile. The main dynamical observed outcome is the mobility 
of the dark component of the soliton. This implies a break-up of the structure; however, 
there are two realizations thereof depending on the symbiotic or non-symbiotic nature of the 
state. In the case of a non-symbiotic solitons (i.e., for NSBBS and for NSDBS), the bright 
component remains at rest forming a genuine bright gap soliton. On the other hand, in the 
symbiotic solitons, this is impossible due to the non-existence of a bright waveform of this 
type. Hence, most of the bright component energy moves towards the opposite direction of 
the dark component in the case of the bubble (i.e., for SEES) while part of the energy moves 
with the dark component. For the SDES, most of the energy appears to move together with 
the dark component. A summary of this scenario is shown in the panels of Fig. [8| To 
indicate the growth rates and unstable eigenmodes of the solutions dynamically followed in 
Fig. [8| we show in Fig. [9] their respective spectral planes. It is worth remarking that, in 
most cases, the instabilities are of exponential and oscillatory type, except in the case of 
SEES, where most of the instabilities are purely oscillatory. On the other hand, to connect 
these results with the experimental motivation of Section [TT| let us point out that for the 
NSEES considered therein the increase of the power is tantamount to a larger instability 
growth rate and hence the observation of the mobility of the dark component, while the 
bright one forms a genuine gap soliton in agreement with our numerics (top panel of Fig. |8]). 
This repulsive effect between the two components is also evident through the blocking of the 
bright channel and the restoration of the bubble at the center, while the reintroduction of 
the interaction between the beams naturally and reversibly reinstates the repulsive bubble 
mobility effect. 



V. CONCLUSIONS AND FUTURE CHALLENGES 

In the work presented in this paper, we have considered the case of two-component dark- 
bright type solitary wave states in defocusing photorefractive waveguide arrays. Motivated 
by experiments in LiNbOs arrays, which illustrated a bubble-type soliton state in one com- 
ponent coupled to a bright gap solitary wave in the second one, we delved into a theoretical 
examination of the different composite states that can emerge in this system. In particu- 
lar, we revealed the potential for four distinct types of waves, namely non-symbiotic and 
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symbiotic, dark-bright and bubble-bright ones. We numerically revealed (within a model 
benchmarked against the linear band structure) the persistence boundaries of such solutions. 
We also analyzed their linear stablity which exhibits a typically weak instability in all of 
them (with fairly small growth rates). This instability is so weak that it permits, appar- 
ently, the experimental observability of the states. Nevertheless, in suitable regimes even 
the experimental dynamics manifests the potential break-up of the composite states. 

Naturally, this investigation paves the way for numerous additional studies. On the 
one hand, from an experimental viewpoint it would certainly be interesting to identify the 
other proposed structures. On the other hand, investigating interactions of such composite 
structures would also offer relevant insights as was done experimentally e.g. with simpler 
states in j22|, or as was done numerically in [40| and in different (BEG) dark-bright contexts 
in [36| [38] . Generalizations of such states in two-dimensional waveguide arrays with the 
formation of vortex-bright states |39j or of genuinely discrete variants thereof [13] would 
also be an exciting theme for future investigations. 
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FIG. 8: Top row: propagation of the squared modulus of the electric field for the dark (left) and 
bright (right) components of a NSBBS with fid = 0.48 [Aneff,d = 2.13 x 10"^] and fit = 0.635 
[Ancff,b = 2.82 X 10"^]. Second row: propagation of the squared electric field for the dark (left) 
and bright (right) components of a SBBS with fid = 0.48 [Aneg^d = 2.13 x 10~^] and /Lt;, = 0.67 
[Anefr,b = 2.98 x lO^^]. Third row: propagation of the squared modulus of the electric field for the 
dark (left) and bright (right) components of a NSDBS with fi^ = 0.32 [Aneg^d = 1-42 x 10~^] and 
fj-b = 0.64 [Arieff^b = 2.84 x 10~^]. Bottom row: propagation of the squared electric field for the 
dark (left) and bright (right) components of a SDBS with = 0.31 [Arieg^d = 1-38 x 10^'^] and 
= 0.66 [Anefr,b = 2.93 X 10-3]. 
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FIG. 9: Stability panels for the solitons of Fig. [8j namely for the NSBBS in the top left, the SBBS 
in the top right, the NSDBS of the bottom left and SDBS of the bottom right. 
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